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The Galileo low-gain antenna mission has a severely rate-constrained channel 
over which we wish to send large amounts of information. Because of this link 
pressure, compression techniques for image and other data are being selected. The 
compression technique that will be used for images is the integer cosine transform 
(ICT). This article investigates the compression performance of Galileo’s ICT al- 
gorithm as applied to Galileo images taken during the early portion of the mission 
and to images that simulate those expected from the encounter at Jupiter . 


I. Introduction 

The Galileo low-gain antenna (LGA) mission will rely 
heavily on source coding, i.e., data compression, to maxi- 
mize the return over the severely rate-constrained chan- 
nel. A major portion of Galileo’s information return 
from Jupiter will be images from the solid-state imaging 
(SSI) camera. Image frames contain either 800 x 800 or 
400 x 400 8-bit pixels, and would require 5.12 Mbits or 
1.28 Mbits, respectively, to transmit uncompressed. Im- 
age compression gains of 10 dB and higher are obtainable 
through known coding techniques such as the discrete co- 
sine transform (DCT), which introduces hardly noticeable 
losses in data fidelity. However, Galileo cannot directly 
utilize off-the-shelf compression algorithms due to its lim- 
ited onboard computing capability. Instead, Galileo will 
use a computationally tractable algorithm called the in- 
teger cosine transform (ICT). The ICT operates and per- 
forms comparably to the DCT. This article investigates 


the compression performance of Galileo’s ICT algorithm, 
as applied to actual Galileo images and to simulated im- 
ages constructed to mimic the conditions that the Galileo 
camera is expected to encounter at Jupiter. 

The amount of compression achieved by the ICT is ad- 
justable by selecting the degree to which the transform 
coefficients are quantized. Coarser quantization generally 
produces higher compression and larger errors in the re- 
constructed image. Thus, the study of ICT compression 
performance is in fact an analysis of the possible trade- 
offs between compression of a source image and the result- 
ing distortions in the reconstructed image. It is impor- 
tant to understand the extent to which the compression 
is predictable by using anticipated image characteristics, 
and how compression versus distortion varies with the se- 
lectable quantization step size. The overall statistical vari- 
ability of the compression performance is an important 
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factor that influences link usage planning. This article is 
a preliminary attempt to examine these issues. 

The ICT is described in [1] and [2], and the specific 
application of the ICT to the Galileo LGA mission is dis- 
cussed in [3]. For a given quantization step size, quanti- 
zation matrix or template, Huffman table, and, of course, 
image, the ICT produces a compressed image. The com- 
pression ratio indicates the ratio of the compressed image 
size to the original image size. Distortion indicates the dif- 
ference between the reconstructed image and the original 
one. Neither the compression ratio nor the distortion is 
input for the ICT. This article describes the behavior of 
the compression ratio and distortion as a function of the 
quantization step size, q ) and any prior information about 
the image. 

An additional level of control may be achievable by 
selecting one of a few specially constructed quantization 
templates and using a corresponding Huffman table. For 
this article, the algorithm uses only a template for uniform 
quantization and a Huffman table that has been trained 
for the template on space images [3]. The tailoring and use 
of different quantization matrices is the subject of ongoing 
work, and will not be discussed here. 

II. The Images in the Data Set 

The ICT has been tested on various images. Ideally, the 
test images would be ones of Jupiter taken by Galileo’s SSI 
camera. Instead, the available images are ones taken by 
Galileo of Venus, the Earth, the Moon, and black sky, and 
ones taken by Voyager of Jupiter and a few of Jupiter’s 
satellites. There are four image sets provided for this 
study by JPL’s Multimission Image Processing Labora- 
tory (MIPL). 

(1) The Earth 1 images are Galileo images from volumes 
2 through 6 of the Earth 1 CD ROMs. They con- 
sist of 970 raw.cal, 1798 Earth, 81 Venus, and 773 
Moon images for a total of 3622 images. A randomly 
selected one-third of the Earth 1 images has been 
used for this article to identify trends and perfor- 
mance predictors. The other two-thirds of the data 
has been reserved for future determination of the ef- 
fectiveness of the predictors and the consistency of 
the trends. 

(2) The 13 selected Galileo SSI images consist of 12 re- 
cent Galileo images and 1 of the Earth 1 images se- 
lected by MIPL to illustrate a typical range of image 
types. These images are designated by names begin- 
ning with rq534. 


(3) The 19 simulated Galileo images were obtained from 
Voyager images of the Jovian system by removing 
Voyager camera characteristics and introducing the 
Galileo SSI camera signature. These images are des- 
ignated by names r.l, r.2, . . . , r.19. 

(4) The 16 radiation-noise- added simulated Galileo im- 
ages are 7 of the simulated images and 9 of the se- 
lected images with 1 of 4 available noise frames (A, 
B, C, and D) added by MIPL to simulate the radia- 
tion effects at Jupiter. These images are designated 
by names beginning with rq538. 

Some examples of the selected Galileo images are shown 
in Fig. 1. The smaller data sets are useful for illustrating 
different behaviors but do not provide sufficient data for 
analysis. The Earth 1 images provide a larger number of 
images for drawing statistical inferences. 

Images were compressed using different values of </, and 
reconstructed. Quantization step sizes 2 through 18, 20, 
and 24 were tested on all images. The compression ratio, 
root mean squared error (rmse), error variance, and error 
mean were noted along with image statistics such as the 
mean, minimum and maximum pixel values, the sample 
zeroth-order entropy, and the sample difference entropy. 
From this information and some a priori information about 
the image, patterns in the compression and distortion per- 
formance have been sought. 


III. Performance as a Function of 
Quantization 

This section shows the relationships between compres- 
sion ratio and image distortion for the ICT as the quanti- 
zation step size is varied. 

A. Compression Behavior 

The compression ratio achieved by the ICT depends on 
the image. This section shows the gross behavior with only 
the small representative image sets. Scatter plots illustrate 
the compression behavior as a function of the quantization 
step size, the spread of compression ratios for different 
images, the limiting compression for very large and the 
effect of the simulated radiation noise on the compression. 

For the Huffman codes used in this work, there is a 
limit of 102.4 on the compression ratio achievable by the 
8x8 ICT. This is because an all-zero quantized block 
in the transform domain can be coded in 5 bits, and an 
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uncompressed block is 8 x 64 bits. In order for a quantized 
block to be all zero, either the original block had to be 
nearly black or the quantization step size had to be very 
large. In either case, the reconstructed block would be 
entirely black. Among the images tested, the nearly black 
ones approached this limit with the quantization step sizes 
tested. 

Figure 2 has three scatter plots showing the compres- 
sion ratio versus quantization step size for the images in 
each of the three small data sets. The plots in Fig. 2 
show that the compression trend is typically monotonic, 
increasing with q , but that it starts to level off for higher 
q. In most instances, as q increases, the spread on the log 
scale of the compression ratio decreases; this is consistent 
with the upper limit on compression ratio. The unlabeled 
scatter plots also show that without taking into account 
information about the specific image, it is difficult to pre- 
dict the compression. 

To illustrate the effects on compression of adding the 
simulated radiation noise to images, Fig. 3 has one plot for 
each of the four noise frames used. The plots show how 
the noise changed the compression ratio for the individual 
images. Solid symbols are used for images in the sets of 
selected Galileo SSI images or simulated Galileo images, 
and hollow symbols with the same shape are used for the 
corresponding images with radiation noise added. These 
plots seem to indicate that the noise frames dominate the 
compression performance; this is more apparent for frames 
A and B and less so for D. 


B. Distortion 

The rmse is a measure of the pixel-by-pixel difference 
between the original image and the reconstructed image. 
It is not the best measure of the distortion of an image, 
but it is easily quantifiable and commonly used as a guide- 
line. A better measure of the distortion is how it affects 
the conclusions and analysis of the images sent back from 
Jupiter, but this is hard to quantify. The SSI team, con- 
sisting of scientists who will be using Galileo’s images, has 
been meeting with both JPL and the Ames Research Cen- 
ter to try to determine what types of distortion are of more 
concern than others. The appropriate distortion measure 
for this analysis is not the subject of this article, and for 
now the rmse will serve as a rough guideline. 

For a given image, the rmse typically increases as the 
quantization step size increases. The three scatter plots 
in Fig. 4 show rmse versus quantization step size for the 
three smaller image sets. Although there is a monotonic 


trend, there are examples of images for which the rmse 
is unexpectedly high for a particular q or q* s. These high 
rmse situations may be avoided most of the time by careful 
selection of the quantization step size. 

C. Selection of Quantization Step Size 

For a particular image, there may be a quantization 
step size that is inferior, meaning that for that image 
there is another quantization step size that gives both bet- 
ter compression and better rmse. An inferior q causes 
a non-monotonic spike on the compression versus distor- 
tion curve. The plot in Fig. 5 shows the distortion ver- 
sus compression ratio for the selected Galileo SSI im- 
ages in Fig. 1, namely rq534.1itn, rq534.ausvn, rq534.1un, 
rq534.pct, rq534.ear9n, and rq534.gas. The solid lines 
show the curves traced out by the rmse and compres- 
sion ratios at all of the tested q’ s; the dashed lines show 
how the curves would be smoother and more predictable 
if the inferior q y s could be predicted and avoided. Images 
that are predominately black sky, such as the image of 
Gaspra, rq534.gas, exhibit the most dramatic rmse fluctu- 
ations with q. 

In order to get a reasonable profile for quantization step 
sizes that are likely to be inferior, the larger image set is 
necessary. The number of images for which a particular q 
was inferior was tallied over the randomly selected third of 
the Earth 1 images. Figure 6 is a histogram showing the 
frequency that different q's were inferior for the raw .cal or 
black sky images, and for the non-raw_cal images that are 
of the Earth, Moon, and Venus. The histogram shows that 
some g’s are inferior more often than others, but sometimes 
only for certain subsets of the images. For instance, q = 6 
is inferior for a quarter of the raw_cal images, but not 
often for the interesting ones. For q = 13, we see the 
opposite trend. Further study is important and should 
include a measure of how bad a q is for a particular image 
(as opposed to simply counting how often q occurs). Also, 
if the images will be classified into types that give improved 
prediction of the compression, then the study of which q } s 
are inferior should be redone for each subset. 


IV. Performance as a Function of Image 
Characteristics 

Using the randomly selected one-third of the Earth 1 
images, correlations were done between the camera set- 
tings and the compression ratio, the compressed image 
size, and the rmse in order to direct the study towards 
likely parameters for predictions, or at least away from 
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parameters that appear to be uncorrelated. The camera 
settings were not continuously valued, and so were not 
ideal for a correlation analysis. A portion of the correla- 
tion matrix of interest is shown in Fig. 7. 

The filter width has a small amount of correlation with 
the performance. The middle wavelength of the camera fil- 
ters seemed independent of the compression and the rmse. 
By dividing the selected data set into subsets defined by 
the different camera settings, statistics for each subset 
could be found. Scatter plots of compression versus rmse 
show the trends for different subsets. However, an analysis 
of variance test with a level of significance of a =: 0.05 in- 
dicated that the subsets probably come from distributions 
with the same mean. 

The black sky images may need to be treated separately 
because compression performance effects that are second 
order for more interesting images are apparently first order 
for these images. For instance, the raw.cal images that are 
mostly black sky tend to cluster according to the gain state 
on the camera (not shown in this article). The effect is less 
noticeable for images with less black sky. 


V. Overall Statistical Variability 

With so many images in the Earth 1 image set, the 
compression and rmse distortion behavior can be displayed 
in histograms where the range of values is divided into 
bins and in cumulative distribution type plots. Figures 8 
and 9 show the histograms and cumulative distributions 
of both compression ratio and rmse for the selected third 
of the images using q — 8. On the cumulative plots there 
are three curves: one is for the Earth, Moon, and Venus 
images; one is for the raw .cal images; and one is for all of 
the images. There are separate histograms corresponding 
to each curve on the cumulative plots. Histograms of this 
sort for different quantization step sizes can assist Galileo 
mission planners in deciding how many pictures to attempt 
to send back, and at what risk this can be done. 


VI. Conclusions 

This preliminary study has generally validated the as- 
sumptions regarding high image compression gains that 
the Galileo LG A mission is counting on. But it has also 
flagged several areas of concern that require further study. 

The scatter plots of Figs. 2-4 and the histograms of 
Fig. 8 show that compression ratios of 10:1 to 40:1 are 
not unreasonable to expect at moderately small values of 
rmse. On the other hand, there is a statistical variability 
of the achievable compression ratio on the order of 5 to 
10 dB from image to image, and slightly less variability 
in the rmse. The variability of compressibility can greatly 
complicate the job of planning how many images can be 
sent back and how large buffers have to be in the onboard 
computer. 

The study of correlations between the known camera 
settings and the resulting compression performance sug- 
gested no significant correlation that could be used to ap- 
preciably improve the predictability of performance. 

There needs to be more analysis of what can be cal- 
culated or predicted about each image to aid downlink 
planning and, perhaps, compression control. For instance, 
better predictors might be obtained from anticipatable or 
measurable image characteristics, such as the percentage 
of black sky expected to fill the image frame. Alternatively, 
an onboard calculation could be used to adjust the quanti- 
zation step size over small portions of the image in order to 
target a particular compression ratio. An example of such 
a calculation of local image statistics is an estimate of the 
entropy using the buffered portion of the image awaiting 
compression. Such ideas are under examination. 

Further improvement in the reconstructed images can 
be achieved by developing more sophisticated reconstruc- 
tion techniques to exploit continuities across block bound- 
aries and other spatial correlations in the image. There 
is ongoing work to develop special purpose quantization 
matrices for improved image quality at good compression 
ratios. 
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Fig. 2. Compression behavior for the small image sets: (a) compression ratio for selected Galileo SSI images, 
(b) compression ratio for images with radiation noise added, and (c) compression of simulated Galileo images 
from Voyager Images. 
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Fig. 4. RMSE distortion for the small image sets: (a) selected Galileo SSI images, (b) images with radiation noise 
added, and (c) simulated Galileo images from Voyager Images. 
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Fig. 5. Effects of inferior quantization step sizes on distor- 
tion versus compression curves. 




Fig. 6. Frequency of occurrence of inferior quantization step sizes over randomly selected images from the Earth 1 image set 318 raw_cal images and 

931 Earth, Moon, and Venus images were used to generate this histogram. 
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Fig. 7. A portion of the correlation matrix. 



FREQUENCY OF OCCURRENCE FREQUENCY OF OCCURRENCE FREQUENCY OF OCCURRENCE 








FRACTION OF IMAGES WITH WORSE COMPRESSION 





